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We show that for a general system of N s-wave point scatterers, there are always N eigenmodes. 
These eigenmodes or eigenchannels play the same role as spherical harmonics for a spherically 
symmetric target — they give a phase shift only. In other words, the T matrix of the system is of 
rank N and the eigenmodes are eigenvectors corresponding to non-0 eigenvalues of the T matrix. 
The eigenmode expansion approach can give insight to the total scattering cross section; the position, 
width, and superradiant or subradiant nature of resonance peaks; the unsymmetric Fano lineshape of 
sharp proximity resonance peaks based on the high energy tail of a broad band; and other properties. 
Off-resonant eigenmodes for identical proximate scatterers are approximately angular momentum 
eigenstates. 

PACS numbers: 03.65.Nk, 34.10.+X, 11.55.-m 



I. INTRODUCTION 

Scattering of waves from a group of scatterers is rich 
in interesting physical phenomena. Examples are elec- 
trons or phonons scattering from defects or impurities in 
crystals and light scattering from conjugated molecules 
or biological complexes. 

Multiple scattering effects are nontrivial, especially 
when scatterers ar e pla ced close, i.e., inside each other's 
effective radius s/cr/w, where a is the total cross sec- 
tion of a single scatterer. If the scatterers are resonant, 
their effective radius can be much larger than their phys- 
ical size or force range. A classical example is the case 
of weak, fixed frequency sound incident on two proxi- 
mate, small identical air bubbles in water 0,0,0). Quan- 
tum mechanically, for two scatterers placed together well 
within the resonant wave length of scattered particle, ex- 
tremely narrow proximity resonance can appear [Q, 
Three or more proximate scatterers lead to related ef- 
fects. 

In this paper we have formulated a novel general 
method, eigenmode expasion, for quantum scattering of a 
system of s-wave point scatterers of any number and geo- 
metric configuration. The eigenmode expansion is based 
on Green function and is equivalent to solving Lippmann- 
Schwinger equation exactly — taking all orders of multiple 
scattering into account. 

It is obvious from eigenmode expansion approach that 
the T matrix of such a system is of rank N and we are 
able to construct the T matrix explicitly. Many features 
of the spectrum (total scattering cross section as a func- 
tion of energy) of identical proximate scatterers can be 



explained simply by this approach. For example, we can 
naturally explain the position, width, and superradiant 
or subradiant nature of resonance peaks, and the un- 
symmetric Fano lineshape of sharp proximity resonance 
peaks based on the high energy tail of a broad band. 
We are also surprised to find off-resonant eigenmodes of 
randomly distributed identical proximate scatterers ex- 
hibit great symmetry and regularity. Similar systems of 
N interacting resonances have been previously studied, 
e.g. H [jj. Though our eigenmode expansion approach 
is of course equivalent to standard scattering theory, it 
gives new physical insights, as we demonstrate below. 

In addition, the quantum scattering of N proximate 
s-wave scatterers is related to the superradiance prob- 
lem H [|. Again, our eigenmode expansion gives new 
insights to superradiance. 

II. REVIEW OF QUANTUM MULTIPLE 
SCATTERING 

We start by introducing our model for s-wave point 
scatterers. Consider a single elastic scatterer placed at 
r = 0, The asymptotic form of the total wave function is: 

vb % {f) ^ e i% - f + f {%,?)- — for r -> oo. (1) 

The first term on the RHS is the incoming plane wave, 
and the second term is the scattered wave. The angular 
dependent scattering amplitude f{k,f) satisfies the opti- 
cal theorem: 4-7r/fc Im/(fc, f = k) = <r(k), where the total 
cross section a(k) = J |/(fc,f)| 2 dr. The T matrix of the 
scatterer in the subspace of energy shell \k\ — k, Tfe, can 
be defined as an operator acting on functions of f: 
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We assume the force range of the scatterer is small com- 
pared to the wavelength of the scattered particle, so 
it only scatters the s-wave component (jo(kr)) of the 
incoming wave. The now angular independent scat- 
tering amplitude is related to the phase angle S(k) by 
f(k) = ±sin<5(fc)e ll5(fc) . The T matrix is of rank 1. T k 
has only one non-0 eigenvalue sin S (k)e i5 ^ , whose corre- 
sponding eigenvector is the Oth order sperical harmonic. 

A simple scatterer with one internal Brcit-Wigner type 
resonance can be modeled as: 

f{k) = -\—^ ■ (2) 

The total cross section has a Lorentzian shape (dashed 
line in Fig. [l]a): 



u ^> k 2 {E~E,) 2 +^ 

where E = k 2 /2 and E = fcg/2. Resonance corresponds 
to the pole of f(E) on complex E plane, at Eq — 270, 
where the "scattered wave" can exist without the incom- 
ing wave. Physically Eq and 70 correspond to position 
and half width of the resonance peak. 

We further assume the force range of the scatterer is 
so small that it can be considered as existing only at 
one point. In this point scatterer model, Eq. (Q) holds 
not only asymtotically, but over all space. Since any 
incoming wave 4>{r) of energy E = k 2 /2 can be written 
as a superposition of plane waves with |fc| = k, the total 
wave function can be written as: 

^(f) = 0(f)+^(O)/(fc)G fc (r), (3) 

where the free space Green function Gf.{r) = e lkr jr. The 
only significant difference is the extra factor <j)(0), the 
amplitude of the incoming wave at the point scatterer. 

We are now ready to tackle the multiple-scatterer prob- 
lem . The total wave function for a system of N scat- 
terers fixed at positions fi, rjy is: 

N 

1>(f) = + 1>i (rDfiGi (r) , (4) 

where /$ is the scattering amplitude of the i-th scatterer, 
Gi(r) is the free space Green function from Gi(f) = 
e lk \ r ~ Ti \ j\r — ri\. and for simplicity we have omitted all 
implicit k (energy) dependence. The 0(0) in Eq. (J3J) is 
replaced by the amplitude at the i-th scatterer of the i-th 
effective incoming wave ipi(r) — the sum of the incoming 
wave and waves scattered by all scatterers except by the 
i-th scatterer itself: 

Ur) = m + E MrDfjGjtf. (5) 

The above equation provides a system of linear equa- 
tions for ipi(fi). We can express them explicitly in matrix 



form. Define vector for incoming wave <j> and vector for 
effective incoming waves ip as: 4n = (f>(fi), ipi = tpi(fi) 
(i = l, — ,N). 
Define the NxN free space Green matrix G as: 

r = / G i(fi) for i ¥"3, 
Uy - \ for i = j, 

and matrix F as: 

F = diag{fx, / 2) • • • , /jv}- 
Substitute r with fi in Eq. (||), we get 

$ = Mip, or %p = MT X 4>, (6) 

where 

M = 1 - GF. 

The physical implication of Eq. (^) is clear. Formally 

M" 1 = 1 + GF + (GF) 2 + ---, 

each term (GF)' (I = 0, 1, • • •) means the incoming wave 
is scattered I times, since F means scattering once by one 
scatterer, and G means free propagation from one scat- 
terer to another. An integral equation similar to Eq. (^) 
is known in general scattering theory as the Lippmann- 
Schwinger equation [JlTJ . Resonances occur at poles of 
M -1 (i?) on complex E plane. 

The scattered wave sc (r) depends only on <fi, the am- 
plitudes of the incoming wave at scatterers. Asymptoti- 
cally 

4>sc = (J2 tife-™*)— for r - 00. (7) 

i=l 

If the incoming wave is a plane wave e tk ' r , the sum in 
parenthesis of Eq. (^) is nothing else but the scattering 
amplitude f(k,f). 

III. FORMULATION OF EIGENMODE 
EXPANSION 

For a single scatterer, for incoming wave jo(kr) (= 
(G(r) — G*(r))/2ik), the scattered wave is proportional 
to G(r) — the outward part of the incoming wave — with 
a phase factor. Is there any incoming wave for multiple 
scatterers such that the scattered wave is proportional 
to the outward part of the incoming wave with a phase 
factor? This is equivalent to looking for eigenvectors of 
the T matrix of the whole system. In the spherically 
symmetric case, diagonalization of T matrix is essentially 
spherical harmonic partial wave expansion. Some text- 
books have discussed the formal diagonalization of the 
T (or S) matrix for a general system, and its applica- 
tion to scattering of particles with spin and scattering 
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reactions Q]. We have found, as will be shown immedi- 
ately, that for a system of s-wave point scatterers of any 
number and geometric configuration, the T matrix can 
be analytically represented under some special basis, and 
can be diagonalized easily. 

From Eq. (§), we see the scattered wave of multiple 
scatterers is a superposition of Gi(r)'s, so we try to write 
the incoming wave as a superposition of jo(k\r — r^|)'s: 

m = EfeioWr- n\) =Egtf g<( °~ Jfe g?(f) , (8) 

i=l i=l 

or = JgJ, where q~j — (qij, • • • , qNj) T , and the N x N 
spherical Bessel matrix J is defined as: 

Jij =jo{k\fi -fj\). 

We have ip = M _1 </>, and from Eq. (Eh, we get 



JV 



ac (f) = 2 Gi(f) [PM" 1 J^] i = ~ £ a(r)[TJ^] 

»=1 i=l 

where matrix T is defined as 

T = fcFM" 1 . 
If q~j is an eigenvector of TJ: 

TJq} = \jq~j (j =1,---,N), 
the scattered wave 

N 



X 



( i>j 



is proportional to the outward part of the incoming wave. 
The angular dependence of the incoming and scattered 
wave is 



N 



Qijt 



-ikf-Ti 



(9) 



Each q~j or yj(f) can be regarded as an eigenmode or 
eigenchannel of the system, is the relative amplitude 
of the scattered wave from the i-th scatterer. Eigenvalue 
Xj is related to the phase angle of the j-th eigenmode Aj 
by 

Xj = sin Aje 1 ^. 

Eigenmodes have some basic properties that can be 
proved straight forwardly. (1) q~j can always be chosen 
real due to time reversal symmetry: 

Qij = Qij- 

(2) The eigenmodes obey orthonormal relations upon 
normalization: 

qfjq'j = J^kj, or J y*{f) yj (f)df = 



(3) Aj's and 2/j(r)'s are actually eigenvalues and eigen- 
vectors of the T matrix: 

T Vj (r)= sin Aje iA i yj (f). 

It is clear now that yj(f)'s here play the same role 
as spherical harmonics in spherically symmetric case. 
Spherical harmonics are eigenchanncls under special con- 
ditions, namely, when the system is spherically symmet- 
ric. 

Any incoming wave can be written as: </>(r) = 
J2iLi c ijo(k\r — + 4>'{r), where 4>'(r) is orthogonal 
to each jo^l^ - Immediately we see the 4>'(f) part 
does not scatter at all, because each scatterer only scat- 
ters s-wave and <j)'(r) is orthogonal to each jo(k\r — 
For a single scatterer, the T matrix is of rank 1. For a 
system of N scatterers, the T matrix is of rank N, since 
the system only scatters incoming waves belonging to the 
iV-dimensional function space the sum in Eq. (||) spans. 

There is a simpler explanation for the rank of the T ma- 
trix. Recall from section || that the scattered wave 4>sc('r) 
depends only on 4>, amplitudes of the incoming wave 
at scatterers. For all the infinitely degenerate incom- 
ing waves of energy E — k 2 /2, we can have a basis in 
which all but ./V waves have amplitude at all fj's. Con- 
sequently, the system only scatters the N waves which 
have non-0 amplitudes at f^'s, and does not scatter the 
rest at all. 

In eigenmodes yj{r) (plus all other linearly 
independent functions of f orthogonal to all 
j/j(f)'s) basis, the T matrix has a diagonal form: 

T = ^a 5 {sinAie iA V--,sinAjve lAjv ,0,---}. The 
eigenmodes are actually eigenvectors corresponding to 
non-0 eigenvalues of the T matrix. In fact, matrix T 
is representation of T matrix in the (nonorthogonal) 
bj(f) = —^J2iLiJij e ~ lkr Ti basis (or JTJ is represen- 
tation of T matrix in the —h==e~ ikf ' r * basis). This can 

in 

be verified if you view Eq. (0) as transformation of basis. 

We can use eigenmode expansion to calculate the to- 
tal cross section for an incoming plane wave. Optical 
theorem can also be verified. The total scattering cross 
section not only depends on the energy (or k), but also 
the direction k of the incoming wave . This is what we 
would expect for a system without spherical symmetry. 
The average of total cross section over all possible in- 
coming directions for a given energy is independent of k. 
From now on, we will use "total cross section" in this 
sense. The total cross section 
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is a sum of contributions from each eigenmode Q. 

Numerical evidence shows that for scatterers with the 
properties of Eq. (^), each phase angle Aj monotonically 
increases from to n as the energy goes from — oo to 
oo. Therefore each eigenmode will contribute a resonance 
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peak to total cross section (if the peaks do not overlap). 
The position of the peak is at where Aj = ir/2, and the 
width is determined by how fast Aj passes through w/2 
and nearby region as a function of energy. For a system 
of N scatterers, there will be N resonance peaks in total 
cross section, whose positions and widths are consistent 
with the N poles of M.^ 1 (E) on complex E plane. 



IV. APPLICATIONS OF EIGENMODE 
EXPANSION TO PROXIMITY RESONANCE 

First we consider the simplest case of multiple scatter- 
ing: two identical scatterers separated by distance d. It 
is easy to verify that F, M, and J all have the form of 
a b 
b a 

ing f(k) explicitly, we can immediately say the two (un- 
normalized) eigenmodes are: 



(a) 



so does TJ 



fcFM _1 J. 



Thus without know- 



Oi 



and q2 



1 

-1 



(10) 



corresponding to the two scatterers oscillating perfectly 
in and out of phase with the same strength. Eq. ( |To| ) is 
valid at any energy. Fig. [l] shows an example of the total 
cross section and phase angles as a function of energy. For 
kd < 1, i.e., the scatterers are placed within each other's 
effective radius, the "s-like" symmetric qi mode becomes 
a broad band, and the "p-like" antisymmetric (fit mode 
becomes a sharp proximity resonance peak, both shifted 
from the peak of single scatterer, as shown in Fig. |l|a [Q. 
In the limit of d — ► 0, qi mode rigorously becomes an 
s-wave and qi mode a p-wave. 

The three scatterer case turns out to have a richer phe- 
nomenology. We consider three identical scatterers po- 
sitioned on a straight line separated by equal distance: 
f x = (-d,0,0),f 2 = (0,0,0),r 3 = (d,0,0). Compared 
to the two scatterer case, the difference is that the rel- 
ative amplitudes in each qj are now energy dependent. 
Nonetheless, upon reflection about yz plane, there are 
always two symmetric modes qi, q^, and one antisym- 
metric mode (ft. Fig. [2Jd shows the phase angles as a 
function of energy. They determine the positions and 
widths of the peaks in total cross section (Fig. |]a). The 
eigenmodes at q{ mode's peak (E = 99.182) are: 



Oi 



0.089 
0.106 
0.089 



0,1 = 




03 = 




We see that qi is s-like, q% p-like, and qi, ci-like. Indeed 
numerical evidence shows that at low energies (E < Eq) 
q[, q~2, and q^ converge to s-,p-, and d- waves when d — ► 0. 
The q{ mode gives most of the broad band and the other 
two modes give two sharp proximity resonance peaks. 

If we disregard the continuity of Aj(E), it is deter- 
mined up to an integral multiple of n. For reasons that 
will become clear in a moment, we have drawn Fig. ^|b 
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FIG. 1: Solid lines: (a) total cross section a and (b) phase an- 
gles Aj of each eigenmode as a function of energy of a system 
of two identical scatterers separated by distance d — 0.2/fco. 
Each scatterer is specified by Eq. (Q) with Eo = 100(a.u.) and 
70 = 0.1. For comparison, the total cross section and phase 
angle of a single scatterer are also shown (dashed lines) . 



with restriction -tt/2 < Aj < it/2. At E ~ 100.578, A : 
and A3 undergo an avoided crossing: instead of one going 
sharply from some —a to 7r — a and the other remain- 
ing at 0, Ai goes sharply from —a to 0~ and A3 goes 
sharply from + to 7r — a. This is really an avoided cross- 
ing of eigenvalues Ai and A3 at A ~ 0. Tracking qi and 
93 around the avoided crossing region, we find qi is s-like 
and 93 is rf-like before the crossing (E < 100.578), yet 
qi becomes d-like and O3 becomes s-like after the cross- 
ing (E > 100.578). In contrast, A 2 rigorously crosses 
Ai and A3, since it has a different symmetry. Indeed if 
we remove the symmetry of the system, there is always 
avoided crossing whenever two phase angles tend to cross 
at (Fig. 0b). The 93 mode peak exhibits a Fano line- 
shape (FigT@a) @. The dip around E - 100.578 before 
the peak comes from the avoided crossing at 0. Sin 2 A 
will exhibit a Fano lineshape if A goes from some —a 
to 7r — a, with minimum at A = and the peak still 
at A = tt/2. Therefore, there is always a dip when two 
phase angles tend to cross at 0. The deepness of the dip 
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FIG. 2: (a) Total cross section a and (b) phase angles Aj 
of each eigenmode as a function of energy of a system of 
three identical scatterers placed on a straight line separated by 
equal distance d = 0.2/feo, at (— d, 0, 0), (0, 0, 0), and [d, 0, 0). 
Each scatterer is the same as used in Fig. |l|. In (b), solid 
lines: Ai, dotted lines: A2, dashed lines: A3. 



depends on how narrowly they avoid crossing. 

Numerical evidence shows that for a general system of 
N identical scatterers confined well within each other's 
effective radius, there is always a broad band with res- 
onance energy < Eq and some sharp proximity reso- 
nance peaks based on the high energy tail of the broad 
band (Fig. 0) jl0| . At peak, the broad band corresponds 
to a mode (say, q{) in which all scatterers oscillate in 
phase, i.e., all qu's have the same sign. The peak is 
broad because the waves scattered from each scatterer 
add up constructively. In this sense this is an s-like mode, 
corresponding to the superradiant state in scattering of 
one photon by a collection of atoms |^| . The remaining 
N — 1 modes give N — 1 sharp proximity resonance peaks 
(if they are separate). They correspond to subradiant 
states. The scattered waves of these modes at their re- 
spective peaks have little s-wave component — waves scat- 
tered from each scatterer tend to add up destructively. 
It has been reported that under certain conditions these 
narrow peak eigenmodes correspond to Anderson local- 
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FIG. 3: (a) Total cross section a and (b) phase angles Aj of 
each eigenmode as a function of energy of a system of three 
identical scatterers placed at (— d, 0,0), (— d/3, — d/3, 0), and 
(d, 0, 0) with d — 0.2/fco. Each scatterer is the same as used 
in Fig. |l[ In (b), solid lines: Ai, dotted lines: A2, dashed 
lines: A3. The difference to Fig. ^ is that the scatterer at 
the center is moved so there is no longer symmetry about the 
reflection of xz and yz plane. 



ization [ fl3| . q{ mode is s-like at low energies. As a con- 
sequence of avoided crossing of phase angles, the s-like 
character will shift to different eigenmodes at different 
energies. Each sharp peak corresponds to a phase angle 
to sharply increasing from some —a (0 < a < tt/2) on 
the broad band (corresponding to s-like mode) to 7r — a — 
back onto the broad band with avoided crossing at 0, ex- 
hibiting a Fano lineshape. We may not see the dip if the 
phase angles do not come close enough in avoided cross- 
ing, but in principle a sharp proximity resonance peak is 
always unsymmetric (Fig. ^a (173 peak only), ||a and [i]a). 

It is surprising to find numerically that, even for a 
typical nonsymmetric system of identical scatterers, at 
far off- resonant energies while still maintaining fcr y <C 1 , 
yj(r)'s are very much like linear combinations of spherical 
harmonics of same I, i.e., they almost have well defined 
angular momenta. yj(f)'s also have lowest Vs possible, so 
they are still s-, p-, d-, ■ ■ ■ like. At resonant energies, the 
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FIG. 4: (a) Total cross section as a function of energy for 
a system of seven identical scatterers randomly placed on a 
plane. Each scatterer is the same as used in Fig. |l| The 
positions of the scatterers are shown in (b) . 



eigenmodes have mixed Us in avoided crossing region. 



ical harmoics for a spherically symmetric target. The 
T matrix of the system is of rank N and can be rep- 
resented analytically. The eigenmodes are eigenvectors 
corresponding to non-0 eigenvalues of the T matrix. The 
position and width of resonance peaks are determined by 
phase angles as a function of energy. Each phase angle 
monotonically increases from to 7r as energy goes from 
— oo to oo if each scatterer has one Breit-Wigner type 
resonance. 

For identical scatterers placed well within one wave- 
length of the incoming wave, the far off-resonant eigen- 
modes are always s-, p-, d-, ■ ■ • like with lowest possible 
angular momenta. At resonant energies where phase an- 
gles have avoided crossings, the eigenmodes have compo- 
nents of different Z's mixed. There is always one broad 
band originating from an s-like mode whose peak is 
shifted lower in energy, and N — 1 sharp proximity res- 
onance peaks based on the high energy tail of the broad 
band. The broad band peak correspond to superradiant 
state in scattering of light, and the sharp peaks corre- 
spond to subradiant states. The unsymmetric Fano line- 
shape of a sharp peak comes from a phase angle sharply 
increasing from some —a on the broad band to n — a 
(back onto the broad band) with avoided crossing at 0. 

In the future, we will investigate the generalization of 
eigenmode expansion to systems of multi-channel point 
scatterers. We will investigate further the dependence of 
the spectrum (total cross section as a function of energy) 
on the spatial distribution of scatterers, and statistical 
properties of the spectra of collections of a large number 
of randomly distributed scatterers. We need to extend 
the theory to non-identical scatterers, and to multipho- 
ton emission of collections of atoms (i.e. true superradi- 
ance). 



V. SUMMARY 

We have shown that for a system of N s-wave point 
scatterers, we can always find N eigenmodes. These 
eigenmodes or eigenchannels play the same role as spher- 
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